arXiv:1503.03898vl [astro-ph.HE] 12 Mar 2015 


Submitted to Astrophysical Journal 

Preprint typeset using style emulateapj v. 05/12/14 


THE MASS AND RADII OF STRONGLY MAGNETIZED NEUTRON STARS 

Farbod Kamiab, Avery E. Broderick, and Niayesh Afshordi 
Perimeter Institute for Theoretical Physics, 31 Caroline Street North, Waterloo, ON N2L 2Y5, Canada 
Department of Physics and Astronomy, University of Waterloo, 200 University Avenue West, Waterloo, ON, N2L 3G1, Canada 

Submitted to Astrophysical Journal 

ABSTRAGT 

It has been clear for some time now that super-critical surface magnetic fields, exceeding 4 x 10^^ 

G, exist on a subset of neutron stars. These magnetars may harbor interior fields many orders of 
magnitude larger, potentially reaching equipartition values. However, the impact of these strong 
fields on stellar structure has been largely ignored, potentially complicating attempts to infer the 
high density nuclear equation of state. Here we assess the effect of these strong magnetic fields on 
the mass-radius relationship of neutron stars. We employ an effective field theory model for the 
nuclear equation of state that includes the impact of hyperons, anomalous magnetic moments, and 
the physics of the crust. We consider two magnetic field geometries, bounding the likely magnitude 
of the impact of magnetic fields: a statistically isotropic, tangled field and a force-free configuration. 

In both cases even equipartition fields have at most a 30% impact on the maximum mass. However, 
the direction of the effect of the magnetic field depends on the geometry employed - force-free fields 
leading to reductions in the maximum neutron star mass and radius while tangled fields increase both 
- challenging the common intuition in the literature on the impact of magnetic fields. 
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1. INTRODUCTION 

Observations and theoretical studies of soft gamma- 
ray repeaters and X-ray pulsars point to the existence 
of neutron stars with very high surface magnetic fields 
{B > 10^'^ G), comprising the so-called magnetars (Dun¬ 
can and Thompson 1992; Paczynski 1992; Thompson and 
Duncan 1995, 1996; Melatos 1999). These surface mag¬ 
netic fields are inferred through the observed slowing of 
the stellar rotation, presumed to be a result of the emis¬ 
sion of energy and angular momentum via large-scale 
magnetic fields at the light cylinder, the point beyond 
which they are unable to continue to rigidly rotate with 
the star. This is expected to spin the star down on a 
timescale « P/P ex B~^P^, where P is the spin period. 
Thus, magnetars are universally observed to have long 
periods, roughly 1 s, and thus correspondingly large light 
cylinders cP/ (27r) « 5 x 10"* km. As a result, the implied 
surface fields necessarily rely on a significant extrapo¬ 
lation, and typically assume a dipolar magnetospheric 
magnetic field geometry, necessarily producing a lower 
limit on the surface field strength, which is itself likely 
to be a lower limit on the interior field strengths. 

As a recent example, Gotthclf et ai (2013) monitored 
the temporal and spectral evolution of a pulsar, originally 
discovered by the NuSTAR X-ray Observatory, and from 
the spin-down measurement, inferred a dipole magnetic 
field strength H = 3 x 10^"* G. Magnetars can also be 
observed in the radio band. Follow-up observations of 
the pulsar PSR J1622-4950, discovered by Levin et al. 
(2010) in a survey of radio pulsars with the Parkes 64 
m telescope, show that the pulsar has the highest in¬ 
ferred surface magnetic field of the known radio pulsars 
(H ~ 3 X 10^"^ G), making it the first magnetar discovered 
via its radio emission. A catalog of 26 currently known 
magnetars was presented recently by Olausen and Kaspi 
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(2014). 

The existence of extremely strong magnetic fields ob¬ 
served in magnetars can be explained by a number of pro¬ 
cesses. Neutron stars with strong magnetic dipole fields 
B ~ 10^'* - 10^® G, can form when conditions for efficient 
helical dynamo action are met during the first few sec¬ 
onds after gravitational collapse (Duncan and Thompson 
1992). In addition to differential rotation, convection 
may play a significant role in amplifying the magnetic 
field (Thompson and Duncan 1993). In the context of 
more sophisticated simulations of field amplification in 
non-rotating stellar cores during the collapse and post¬ 
bounce accretion phases of a supernova, including mag¬ 
netohydrodynamics and neutrino transport, it was found 
that initial magnetic field strengths stronger than 10^° 
G could yield magnetar-like final field strengths (Ober- 
gaulinger et al. 2014). 

As neutron stars with very strong magnetic fields ap¬ 
pear to exist in nature, one is tempted to ask how these 
field strengths affect the structure of these stars. This 
is particularly important as the equilibrium mass and 
radius of neutron stars vary based on the nuclear and 
gravitational physics assumed (for example, see Lattimcr 
and Prakash (2001, 2007) for the effects of various nu¬ 
clear equations of state on the structure of neutron stars, 
and see Kamiab and Afshordi (2011); Pani et al. (2011); 
Alavirad and Weller (2013) for effects of modifying gen¬ 
eral relativity). As a result, observational measurements 
of these masses and radii have the potential to con¬ 
strain theoretical models. In particular, observing neu¬ 
tron stars with very high masses is useful, as each set 
of models (nuclear equation of state and gravitational 
model) predicts a maximum mass beyond which no neu¬ 
tron stars would exist. For example, the detection of a 
1.97 ± 0.04 Mq pulsar by Demorest et al. (2010), or the 
measurement of a 2.01 ± 0.04 Mq pulsar by Antoniadis 
et al. (2013) have been used to significantly constrain 
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the viable nuclear equations of state, as well as potential 
modifications of general relativity (Kamiab and Afshordi 
2011). Therefore, before reaching definite conclusions, it 
is crucial to include all the necessary physics, and in the 
case of magnetars, investigate how strong magnetic fields 
affect their mass-radius (M-R) distribution. 

The presence of a strong magnetic field in a neutron 
star can potentially affect the mass and radius in two 
ways: 

1. Locally, a magnetic held affects the nuclear equa¬ 
tion of state (EOS) thus indirectly affecting the 
equilibrium conhgurations of the star. 

2. The magnetic held can affect the structure of the 
neutron star globally by contributing to the hy¬ 
drostatic support against gravity via its stress, 
and the structure of its spacetime via its energy 
(in other words, through the contribution of the 
Maxwell energy-momentum tensor to the Einstein 
equations). 

Broderick et al. (2000, 2002) studied the effects of a 
strong magnetic held on the nuclear EOS. As this work is 
an extension of their studies, we will describe this magne¬ 
tized nuclear EOS in more detail in Section 2. The EOS 
expresses nuclear pressure in terms of the local magnetic 
held strength and nuclear density. Therefore, some as¬ 
sumption about the structure of the magnetic held is 
necessary. Cardall et al. (2001) studied static neutron 
stars with poloidal magnetic helds and a simple class of 
electric current distributions consistent with the require¬ 
ment of stationarity. Considering the global ehect of the 
magnetic held stress in the Einstein equations, and as¬ 
suming a set of nuclear EOS, mainly the one calculated 
by Broderick et al. (2000, 2002), they found that the 
magnetic held increases noticeably the maximum mass. 
In a recent study, Astashenok et al. (2014) have consid¬ 
ered neutron stars with strong magnetic helds (where the 
held strength is a simple parametric function of baryon 
density only) in the framework of f{R) gravity. Using a 
nuclear EOS similar to the one calculated by Broderick 
et al. (2000, 2002), and assuming the global effect of the 
magnetic held pressure, they found that the strong mag¬ 
netic held can increase considerably the maximal mass 
of the star. They also hnd that for large helds, the M-R 
relation differs considerably from that of general relativ¬ 
ity only for stars with masses close to the maximal one. 
In another recent study, Lopes and Menezes (2014) as¬ 
sume a chaotic magnetic held model introducing a vari¬ 
able magnetic held, which depends on the energy density 
rather than on the baryonic density, and based on this 
calculate the mass radius relationship for neutron stars. 

It is important to stress that in the studies of Lopes and 
Menezes (2014), Astashenok et al. (2014) and Cardall 
et al. (2001), the magnetic held is able to provide hydro¬ 
static support in the global sense mentioned above where 
the magnetic held contributes directly in pressure in the 
Einstein equations [for example, Pb = in the case 

of Astashenok et al. (2014) and Pb = 124,7: in the case 

of Lopes and Menezes (2014)]. Although the local effects 
of the magnetic held on the EOS should be considered in 
calculating the M-R relations, it is not certain whether 
the magnetic held can provide global hydrostatic sup¬ 
port for the star. Using the principle of conservation of 



Figure 1. EOS - pressure (P) versus number density of baryons 
(n;,) for different values of the magnetic field strength (f?) assuming 
no anomalous magnetic moments. Strong magnetic fields soften the 
EOS because of Landau quantization. 

total helicity, Broderick and Narayan (2008) developed a 
variational principle for computing the structure of the 
magnetic held inside a conducting sphere surrounded by 
an insulating vacuum. They show that, for a hxed to¬ 
tal helicity, the minimum energy state corresponds to a 
force-free conhguration, which is generically anisotropic. 
If magnetic held lines rearrange to a force-free conhgu¬ 
ration in neutron stars, then by dehnition, they cannot 
hydrostatically support the star, and thus their impact 
is limited to local effects on EOS. 

The purpose of this work is to study the effect of these 
various assumptions (in particular, the force-free model 
of Broderick and Narayan 2008) on the M-R distribution 
of neutron stars. Eor the effects of the magnetic held on 
the nuclear EOS, we base our work on the calculations 
of Broderick et al. (2000, 2002) which will be briehy de¬ 
scribed in Section 2. In Section 3 we will motivate the 
particular magnetic held prohle we consider, based on a 
simplihed model of neutron star formation and magnetic 
held amplihcation that results in a local magnetic held 
strength hxed to a fraction of the equipartition value. 
By varying the fraction of equipartition we are able to 
study the effect of diherent held strengths on the M-R 
distribution. In Section 4 we will calculate the M-R re¬ 
lation of neutron stars, assuming diherent strengths of a 
statistically isotropic magnetic held, based on diherent 
models for the nuclear EOS. In Section 5 we will study 
the structure of neutron stars with a force-free magnetic 
held prohle (which is anisotropic). Einally, we will con¬ 
clude in Section 6. 

2. MAGNETIZED NUCLEAR EQUATION OF STATE 

In this work, we employ the EOS calculated by Broder¬ 
ick et al. (2000), which includes the ehects of very strong 
magnetic helds in multicomponent, interacting matter. 
This is obtained via a held-theoretical approach in which 
the baryons (neutrons, n, and protons, p) interact via the 
exchange of a-uj-p mesons. They consider two classes 
of models that diher in their high densities behavior. 
In these models, the ehects of the anomalous magnetic 
moments of the nucleons are introduced via a covariant 
coupling of the baryons to the electromagnetic held ten¬ 
sor, presumed to be the result of of an appropriate high- 
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Figure 2. Same as Fig. 1, but now including the effect of anoma¬ 
lous magnetic moments. The magnetic field stiffens the EOS due 
to the polarization of the anomalous moments until complete spin 
polarization occurs {B > 10^® G). From that point on, the effect 
of Landau quantization takes over and softens the EOS. 

energy theory. The magnetic fields affect the EOS in two 
ways. In strong fields {B > 10^^ G), the orbital motion 
of the charged particles is quantized (Landau quantiza¬ 
tion), resulting in a reduction in the electron chemical 
potential and a substantial increase in the proton frac¬ 
tion. This produces a softening of the EOS. However, if 
the anomalous moments are considered, the authors find 
that for very strong fields {B > 10^® G), complete spin 
polarization of the neutrons occurs. This produces an 
overall stiffening of the EOS that overwhelms the soften¬ 
ing induced by Landau quantization. These effects are 
illustrated in Figures 1 and 2. 

Broderick et al. (2002) have also studied the effects of 
hyperons on the magnetic EOS. Their EOS with hyper- 
ons, in the absence of a magnetic field is equivalent to the 
GM3 EOS calculated by Glendcuning and Moszkowski 
(1991) where hyperons are introduced as free baryons, 
interacting only through weak interaction so that the 
system maintains nuclear statistical equilibrium. Typ¬ 
ically, the introduction of the additional degrees of free¬ 
dom associated with the new baryon species produces a 
net softening of the equation of state. 

Since the pressure is now generally a function of both 
the matter density and magnetic field strength, explicit 
construction of a neutron star requires some knowledge 
about the magnetic field profile. 

3. A TANGLED MAGNETIG FIELD 

The primary processes responsible for the amplifica¬ 
tion of the magnetic field are believed to be dynamos 
driven by differential rotation and convection Duncan 
and Thompson (1992); Thompson and Duncan (1993). 
These necessarily produce a tangled magnetic field inside 
the neutron star. Since for simplicity we assume spheri¬ 
cal symmetry throughout this work, we must model the 
dependence of the strength of this tangled magnetic field 
on radius. We do this approximately by making a hand¬ 
ful of simplifying assumptions regarding the outcome of 
the dynamo process. 

Various studies suggest that the dynamos naturally 
saturate at locally equipartition magnetic field strengths 
during the initial formation of the proto-magnetar after 


Figure 3. EOS - pressure (P) versus number density of baryons 
{ni,) for different values of the magnetic field strength (B) assuming 
the existence of magnetic anomalous moments. The solid lines are 
the EOS found by imposing the equipartition condition (Equation 
2) using different values of 0. Smaller values of 0 are associated 
with higher values of magnetic field strength. 

the supernova explosion, while different layers of the re¬ 
maining stellar core condense to form the neutron star 
(Chevalier 2005; Naso et al. 2008; Obcrgaulingcr et al. 
2014). That is, we expect 

es oc eg, (1) 

where es is the local energy density of the magnetic field 
and eg is the local energy density of the gas. Since Pg « 
eg/s at this time, this immediately implies 
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( 2 ) 


where /3 is the standard proportionality factor relating 
the gas and magnetic pressures. 

As the proto-magnetar cools, /3 will be approximately 
conserved. During the cooling process the star will main¬ 
tain hydrostatic equilibrium and hence for each spherical 
shell within the star 


Pn 


GM 


P- 


( 3 ) 


where r is the shell radius, M is the enclosed mass, p is 
the matter density, and G is Newton’s constant. As¬ 
suming that the magnetic field has grown sufficiently 
strong to suppress any further convection, M will remain 
constant and the proto-magnetar will contract homolo¬ 
gously. Thus p will grow roughly proportional to dur¬ 
ing the subsequent proto-magnetar evolution, and thus 
Pg oc r~‘^. Similarly, the flux conservation implies B will 
grow oc r“^, and thus B"^/Sn oc r~‘^ as well. Therefore, 
we expect /3 to remain fixed during the proto-magnetar’s 
formation following the quenching of the dynamos. 

This relation naturally provides a profile for the mag¬ 
netic field which is proportional to the pressure of mat¬ 
ter; given a global value of j3 we may identify a unique 
pressure at each density. Examples of the EOS found 
by imposing Equation (2) on the results of Broderick 
et al. (2000) assuming anomalous magnetic moments are 
shown in Figure 3. 

In addition to the modification of the nuclear equa¬ 
tion of state, the strong interior magnetic fields may pro- 
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duce magnetic stresses directly. However, the nature of 
these stresses depend on the local geometry of the field 
as well as its strength. Assuming that the outcome of 
the dynamo processes is a small-scale, tangled field that 
is weakly uncorrelated with the generating currents, the 
resulting stresses may be completely described by a mag¬ 
netic pressure, which after assuming statistical isotropy 
is given by 


Pb 


1 

3 8^’ 


(4) 


The factor of 1/3 arises from the impact of magnetic 
tension; that it must be present is immediately evident 
from the fact that the electromagnetic stress-energy ten¬ 
sor is traceless. In Equation (2), (3 = 1/3 (giving the 
EOS lying between the cyan and blue curves in Figure 
3), corresponds to a neutron star with equal pressure of 
magnetic field and matter. Smaller values of /3, for which 
the pressure of the magnetic field will be larger than the 
one of matter, are not generically expected in turbulent 
MHD dynamos, as magnetic energy is converted back 
into thermal energy through magnetic reconnection. As 
we see in Figure 3, complete spin polarization of neutrons 
(due to their anomalous magnetic moments) occurs at 
(3 = 0.001. Therefore, it is expected that at values of 
/3 equal to, or larger than 1/3, the anomalous magnetic 
moments will not have a considerable effect on the nu¬ 
clear EOS and consequently in the M-R relations. We 
will see this explicitly in the next section. 

We note in passing that the assumptions made above 
regarding the correlation between the generating currents 
and the magnetic field need not be typically satisfied. An 
extreme alternative are the force-free configurations de¬ 
scribed in Broderick and Narayan (2008), with a corre¬ 
sponding impact on the nature of the magnetic stresses. 
However, we will postpone a detailed discussion of this 
until Section 5. 

Before describing our M-R calculation, it is worth men¬ 
tioning that the magnetized EOS of Broderick et al. 
(2000, 2002) describes the nuclear matter inside the neu¬ 
tron star. For the crust, i.e., densities below the nu¬ 
clear saturation density (baryon rest-mass density of 
Pns ~ 2.7 X 10^"^ g cm^), we use the SLy (Skyrme Lyon) 
EOS calculated by Doiichiii and Haenscl (2001), which 
is based on the effective nuclear interaction SLy of the 
Skyrme type, which is useful in describing the properties 
of very neutron rich matter. 


4. M-R RELATIONS WITH AN ISOTROPIC MAGNETIC 
FIELD 

As choosing a magnetic field configuration given by 
Equation (2), allows us to calculate a unique magne¬ 
tized EOS (as seen in Figure 3), we can now proceed 
to calculating M-R relations for neutron stars. As we 
are assuming spherical symmetry and isotropy, the Ein¬ 
stein equations in hydrostatic equilibrium, simplify to the 
Tolman-Oppenheimer-Volkoff (TOV) equation: 


dP 

dr 




M + 47rr 


3 



2GM- 
c^r . 


-1 


(5) 



Figure 4. The M-R relation for the magnetized EOS (Broderick 
et al. 2000) with no crust physics included. Here, the magnetic field 
provides hydrostatic support for the star by entering globally in the 
Einstein equations as Pb = (l/3)(-B^/87r). The calculation has 
been done for different values of jS in Equation (2). Smaller values 
of /3 are associated with higher values of magnetic field strength 
in the neutron star. = 1/3 corresponds to a neutron star with 
equal pressure of magnetic field and matter. 

where c is the speed of light, and e and P are the total 
energy density and pressure: 

e = eg + eB, ( 6 ) 

P = Pg+PB, (7) 

where eg and Pg (the energy density and pressure of the 
gas), are related by imposing Equation (2) (with a fixed 
choice of /3), on the various magnetized nuclear EOS 
models of Broderick et al. (2000, 2002) (as shown in Fig¬ 
ure 3 for one particular nuclear EOS). Where stated we 
include the physics of the crust, by using the SLy EOS 
of Douchin and Haenscl (2001) below nuclear saturation 
density. The energy density of the magnetic field is given 
by cr = H^/Stt, where again B is fixed by the choice of (3. 
The pressure of the magnetic field is given by Equation 
4 (statistically isotropic force-inducing magnetic field). 
The mass enclosed inside radius R is is related to local 
energy density by 

pR 

M{R) = / e(r) Airr^dr. (8) 

Jo 

We integrate Equation (5) numerically, starting from 
the centre of the star and assuming different values of 
central pressure Pg{r — 0). Figure 4 shows the M-R 
relations resulting from this integration in the absence of 
hyperons and crust physics. Increasing the magnetic field 
strength (decreasing j3) in the absence of hyperons has 
only a small effect on the neutron star maximum mass. 
However, larger magnetic field strengths are associated 
with larger neutron star radii. 

Figure 5 shows the M-R relations in the presence of hy¬ 
perons, again ignoring the crust physics. In the presence 
of strong magnetic fields, all of the hyperons are suscep¬ 
tible to spin polarization. As we discussed in Section 2, 
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Figure 5. Same as Figure 4, but now including the effect of hy- 
perons. 

spin polarization counteracts the effects of Landau quan- 
tization, by stiffening the EOS. Although assuming the 
existence of hyperons gives a smaller maximum mass in 
the absence of a magnetic field (/3 —?► oo in Figure 5), in¬ 
creasing the magnetic field strength (decreasing /3), now 
increases the maximum mass. Assuming equal pressures 
of magnetic field and matter (/3 = 1/3), we get a max¬ 
imum mass which is 26% larger than the one with no 
magnetic field (/3 —>■ oo). 

Figure 6 shows similar results, now incorporating the 
physics of the crust by using the SLy EOS calculated by 
Douchin and Haenscl (2001) at densities below nuclear 
saturation density. In this case, the M-R relation found 
by assuming zero magnetic field strength (^ —>■ oo), is 
similar to the well-known GM3 EOS calculated by Glen- 
denning and Moszkowski (1991). A similar increase in 
neutron star maximum mass with magnetic field strength 
is seen here.^ 

We are also interested in the effects of the anoma¬ 
lous magnetic moments of nucleons on the M-R rela¬ 
tions. Figure 7 shows these results for the magnetized 
nuclear EOS with hyperons and crust physics included. 
As we see, the effect of including the anomalous magnetic 
moments is small. We expected this from applying the 
condition given by Equation (2) on the magnetic EOS of 
Broderick ef al. (2000, 2002) with anomalous moments. 
As can be seen in Figure 3, complete spin polarization 
of the anomalous moments (which is where we would ex¬ 
pect their effects to be considerable) occurs at values of 
P ~ 0.001, in other words when the pressure of the mag¬ 
netic field is a 1000 times larger than pressure of matter, 
which are not expected to arise in dynamos, as discussed 
above. Therefore the strongest magnetic field strength 
we considered was one for which the magnetic pressure 
equals the pressure of matter (/3 = 1/3). The effect of 
including magnetic anomalous moments under this mag¬ 
netic field strength is a 2% change in the neutron star 

^ For another recent study of the effects of strong magnetic fields 
on the population of hyperon stars, see Gomes et al. (2014). 



Figure 6. Same as Figure 5, but now also including the crust 
physics via the SLy EOS (Douchin and Haonsel 2001). The black 
curve, corresponding to negligible magnetic field strength, is equiv¬ 
alent to the GM3 EOS (Glondenning and Moszkowski 1991). 



Figure 7. The M-R relation for the magnetized EOS with hyper¬ 
ons (Broderick et al, 2002) and the crust physics included via the 
SLy EOS (Douchin and Haenscl 2001). Here, the magnetic field 
provides hydrostatic support for the star by entering globally in 
the Einstein equations as Pg = (l/3)(B^/87r). The solid curves 
have been calculated assuming no anomalous magnetic moments 
and the dashed curves have been calculated assuming the exis¬ 
tence of anomalous moments. The calculation has been done for 
different values of 0 in Equation (2). Smaller values of 0 are asso¬ 
ciated with higher values of magnetic field strength in the neutron 
star. 0 = 1/3 corresponds to a neutron star with equal pressure 
of magnetic field and matter. The black curve, corresponding to 
negligible magnetic field strength, is equivalent to the GM3 EOS 
(Glendenning and Moszkowski 1991). 
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Figure 8. Change in maximum mass of neutron stars as a function 
of 1//3 (see Equation (2) for definition of /3) for the mass-radius 
relation shown in Figure 6. 


Figure 9. Change in radius of neutron stars as a function of 1//3 
(see Equation (2) for definition of /3) for the mass-radius relation 
shown in Figure 6. 


maximum mass as seen in Figure 7. 

Finally, to better illustrate the effect of the magnetic 
field strength on the maximum mass and radius of neu¬ 
tron stars. Figures 8 and 9 show the maximum mass and 
inflection point in radius (where the second derivative of 
mass versus radius changes sign) as a function of 1//3, 
for the mass radius relationship shown in Figure 6. We 
specify the radius by the inflection point in the M-R rela¬ 
tion, where mass is nearly independent of radius, giving 
a characteristic value across a wide range of neutron star 
masses. We see that these relative changes are well fit by 


cs 0.1 r\ and cs 0.06 (9) 




7?inf 


for an isotropic tangled magnetic field distribution. 


5. M-R RELATIONS WITH A FORCE-FREE MAGNETIC 
FIELD 

Despite the production of small-scale, turbulent fields 
by the dynamos initially, the final magnetic field con¬ 
figuration remains unclear. After the dynamos quench 
and before the formation of a crust that can balance 
shear stresses the neutron star will reconfigure the mag¬ 
netic field geometry via bulk fluid motions on a timescale 
measured in Alfven crossing times, resulting in a linked, 
nearly force-free geometry dictated by the initial mag¬ 
netic helicity (Braithwaite and Spruit 2004, 2006). This 
is a natural consequence of helicity conservation dur¬ 
ing reconnection, corresponding to the minimum energy 
state at fixed magnetic helicity (Broderick and Narayan 
2008), which is fully dehned by 

V X B = oB , (10) 

from which it is clear that j x B = 0 where j is the current 
density, i.e., the Lorentz force vanishes. In the above, a 
is a constant set by the boundary conditions at the stellar 
surface, setting a scale length for the resulting magnetic 
geometries. 


This gives axisymmetric and anisotropic conhgurations 
for the magnetic field. This configuration in the interior 
of the star, in terms of the vector spherical harmonics 
Yi, and is: 


B-E - 


+ 1 ) J • ^ 

- diji{ar) 


Yi 


drrdiji{ar) 


+adiji{ar)^i} , (11) 


where I is the spherical harmonic degree, ji is the spher¬ 
ical Bessel function, and di is a constant giving the 
strength of each harmonic mode. 

To assess the consequence of the force-free condition 
on the hydrostatic equilibrium of the star, we can start 
by looking at the Newtonian Euler equation and impose 
static conditions: 

=-VPg-pV$-bj X B =-VPg-pV$ = 0, (12) 

where p is the density and $ is the Newtonian potential 
given by the Possion equation: 

= 47rGp. (13) 

Therefore, we have: 

VPg = -pV$, (14) 

which is the Newtonian equation for hydrostatic equilib¬ 
rium with no effect of the magnetic field in it. This means 
that in Newtonian mechanics, the force-free condition 
implies that the magnetic field will play no global role in 
supporting the star, even though it might still change the 
nuclear equation of state via the effects mentioned above 
(e.g., Landau quantization, spin polarization, etc.). 

In general relativity, the force-free condition can be 
written in its covariant form: 


= 0 ^ = 0, (15) 
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where F^i, is the electromagnetic tensor. The electro¬ 
magnetic stress-energy tensor is given by 

= (16) 

The full stress-energy tensor is given by: 

= Tgas^'' + Tem^’', (17) 

from which the Euler equation is obtained by = 0. 

The contribution from Tem^'' is generally, after some 
manipulation and application of Maxwell’s equations, 

VuTem^'' = = j^F^'' = 0, (18) 

47r 

therefore as in the Newtonian limit the direct impact 
of the electromagnetic forces vanish, i.e., V= 

y rp _ Q 

v gas — 

Gravity is now described by the Einstein equation: 

= 8^G(Tg„/" + . (19) 

Unlike the Newtonian case, now the electromagnetic field 
can contribute by sourcing gravity. In practice, it is 
necessary to solve for the full magnetic held conhgu- 
ration and metric simultaneously. However, we might 
imagine that since gravity is a long-range force and the 
structure of the initially highly turbulent electromagnetic 
held should exhibit mostly small-scale structure that we 
may spatially average the electromagnetic conhguration. 
That is to solve: 


= SirGiTgas^’' + (TEMn), (20) 

where (...) is a local spatial average, i.e., an average over 
scales large in comparison to the scale of the electromag¬ 
netic held huctuations and small in comparison to the 
gravitational scales. The stress energy tensor for this 
average held is 


/ 

where 

eB = {B^)/ST:, 

pP = (P"}/87r - (P^Pr}/47r, 

pjf^ = (B^}/8n-(B^Bf,}/47r, 

pjf^ = (B^}/87r-(B<^B^}/47r. 

The anisotropy of this average conhguration may be 
quantihed via an anisotropy parameter. 


(TemH = 

V 


A = 


IP. 


W+pW]/2-p^ 


(r) 


(R2)/8^ 


( 21 ) 


An anisotropy parameter A = — 1 corresponds to a 
purely tangential magnetic held conhguration ((i?’’) = O) 
and A = 2 corresponds to a purely radial held ((B'’) = 
(B)). An isotropic average held will have A = 0. 

Figures 10 and 11 show the anisotropy parameter cal¬ 
culated for the average force-free magnetic held conhgu¬ 
ration given by Equation (11) for two spherical harmonic 


A 



Figure 10. Anisotropy parameter for the average force-free con¬ 
figuration for the I = 1 spherical harmonic degree. Units of radius 
are set so that a = 1 in Equation (11). 


A 



Figure 11. Anisotropy parameter for the average force-free con¬ 
figuration for the I = 10 spherical harmonic degree. Units of radius 
are set so that a = 1 in Equation (11). 


degrees of Z = 1 and I — 10. At the centre of the star, 
there is a region with an isotropic average held (A = 0) 
for Z = 1 and as we go to higher spherical harmonic de¬ 
grees A —>• 0.5 at the centre as can be seen in Figure 
11. As we see in both Figures 10 and 11, A drops to 
-1 towards the outside of the star. This matches our 
expectation of held lines closing and the magnetic held 
becoming tangential in the outside. In the calculations 
below, we will refer to the radius of the isotropic core 
with A = 0 as B* (having chosen the anisotropy param¬ 
eter of the hrst harmonic degree, as the magnetic held 
strength is dominated by the hrst harmonic degree at the 
centre). 

To calculate the mass-radius relationship of the star, 
we will assume spherical symmetry and continue using 
an equipartition magnetic held (Equation (2)). In order 
to write the hydrostatic equilibrium equations coming 
from general relativity for the force-free model, we use 
Equation (18) (in other words the fact that the Euler 
equation is simply given by the energy momentum con¬ 
servation of the gas). As an approximation, we will solve 
the Einstein equation given by Equation (20). It is worth 
clarifying that solving Equation (20) alongside the en¬ 
ergy momentum conservation of the gas uTgas^'' = 0), 
would formally require 

V^(Ti5Mn=0, (22) 

while by using the force-free condition Equation (18) we 
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are only assuming that 


= 0. (23) 

As we will show in the Appendix, solving Equation (22) 
alongside with the Euler equation for the gas, would 
uniquely determine the anisotropy parameter. As spher¬ 
ical symmetry only allows isotropy (A = 0) at the centre 
of the star, solving Equation (22) will result in A —)■ 0 
throughout most of the star. This will change the final 
mass-radius relationship results we obtain from using the 
anisotropy parameter coming directly from the force-free 
configuration (satisfying Equation (23)) at most by 9%. 
This anisotropy parameter based on the results shown in 
Figure 10 will have a core with A = 0 of size i?* and in 
our approximation, A will drop to —1 at i?*. 

In dimensions where G = 1 and c = 1, the Euler equa¬ 
tion coming from 'S/i,Tgas^'' = 0 is given by 


where g is 


dP^ 

dr 


-(Cg + Pg) g, 


M -f 471 (Pg -h P^’'^) r3 
r2(l - 2M/r) 


(24) 


(25) 


The enclosed mass is still given by Equation (8), which 
still includes the magnetic contribution to the energy 
density. 

Based on the definition of A and the assumption of 
spherical symmetry {B^Bg) = (B^B^), from which 

^ 1^ ( 26 ) 

Using this and the equipartition assumption the total 
radial pressure becomes 

Pg+P« = (l+i^)Pg, (27) 

requiring only the specification of Pg and A. The pres¬ 
sure and energy density of the gas are related through 
the nuclear equation of state which is itself affected by 
the magnetic field. The resulting mass-radius relation¬ 
ships from this numerical integration are shown in Figure 
12. The solid curves correspond to the force-free mass- 
radius calculation based on the magnetized EOS with 
hyperons (Broderick et al. 2002) and the crust physics 
included via the SLy EOS (Douchin and Hacnscl 2001). 
This calculation has been done for two different values of 
/3. The range in each curve is due to assuming a range 
in P* (radius of the isotropic core with A = 0) from 0 
to 3 km. The previously found mass-radius relationships 
based on an isotropic magnetic field are also plotted for 
comparison (dashed curves). 

In stark contrast to the case of tangled, isotropic mag¬ 
netic helds, a force-free configuration produces a substan¬ 
tial decrease in the maximum neutron star mass. This is 
not surprising - force-free magnetic field geometries con¬ 
tribute to the energy density, and therefore gravity, but 
not to the supporting pressure. This conclusion is only 
weakly dependent on the details of the force-free config¬ 
uration under consideration, showing little sensitivity to 
the size of our assumed isotropic core, P*. However, even 



Figure 12. The M-R relations for the magnetized EOS with hy¬ 
perons (Broderick et al. 2002) and the crust physics included via 
the SLy EOS (Douchin and Hacnscl 2001). The dashed curves are 
for an isotropic magnetic field and the the solid curves correspond 
to the force-free model. Results are shown for two values of (3. The 
range in the force-free (solid) curves is due to assuming a range in 
R* (radius of the isotropic core with A = 0) from 0 to 3 km. 


o 10 
o 

X 

S 

< 


Figure 13. Change in maximum mass (absolute value for the 
logarithmic scale as the mass is decreasing) of neutron stars as a 
function of 1//3 (see Equation (2) for definition of /3) for the force- 
free mass-radius relations shown in Figure 12. 

for large values of /3, e.g., (3 = 1, the magnitude of the ef¬ 
fect is constrained to 20%, similar in size if not direction 
to that associated with the isotropic tangled geometry. 

The dependence on magnetic field strength of the size 
of the effects on mass and radius are illustrated in Figures 
13 and 14. Note that unlike Figures 8 and 9, the 
change in mass and radii are now negative. As before, 
the change in mass and radii are well fit by power laws 
over at least an order of magnitude in /?, for maximum 
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Figure 14. Change in radius of neutron stars (absolute value for 
the logarithmic scale as the radius is decreasing) as a function of 
1//5 (see Equation (2) for definition of f3) for the force-free mass- 
radius relations shown in Figure 12. 


mass and inflection point radius 

~ -0.13 /3-0-9, and ~ -0.03 


M„ 


.^inf 


( 28 ) 

assuming a force-free tangled magnetic field distribution. 


6. CONCLUSIONS 

In this paper, we have demonstrated that the assumed 
geometry of interior neutron star magnetic fields is at 
least as important as the held strengths themselves when 
assessing the impact on masses and radii, capable of qual¬ 
itatively modifying their effect. Highly tangled, isotropic 
magnetic conhgurations do result in an increase in the 
maximum mass, typical of results often found in the 
literature (e.g., see Astashenok et al. 2014). However, 
force-free conhgurations, motivated by dynamical studies 
of magnetic held evolution in proto-neutron stars, pro¬ 
duce a deerease in the maximum neutron star mass. For 
both conhgurations of the equipartition strength helds, 
the magnitude of the impact on the maximum mass is 
limited to 30%. The impact of anomalous magnetic mo¬ 
ments for equipartition magnetic held strengths at best 
modify the maximum mass by 2%, suggesting that they 
are sub-dominant to the global magnetic held geometry. 

This immediately challenges the assertion that mag¬ 
netic helds provide a means to reconcile the recent obser¬ 
vations of very massive neutron stars with unmagnetized 
equations of state whose maximum mass is otherwise pre¬ 
cluded. In the absence of some credible mechanism for 
producing super-equipartition magnetic helds, the maxi¬ 
mum deviation in the maximum mass of order 20%-30% 
means that at best these can impact only marginal cases. 
For example, this is unable to bring the GM3 EOS we 
employed into consistency with the existing 2 Mq mass 
measurements. We expect this to be generic for models 
of neutron star formation in which the interior magnetic 
helds are produced by turbulent processes, providing an 


optimistic upper limit on the mass-limit enhancement 
achievable via the introduction of electromagnetic sup¬ 
port. Of course, the existence of well-motivated magnetic 
held conhgurations in which the maximum neutron star 
mass actually decreases (i.e., force-free conhgurations) 
calls into question the value of invoking magnetic sup¬ 
port altogether in order to increase maximum neutron 
star masses. 

We note that the observed anomalous braking in¬ 
dices for the spin-down of pulsars can be interpreted 
as evidence for sub-surface toroidal magnetic helds of 
10^^ — 10^® G, substantially stronger than the observed 
dipole helds, which are slowly diffusing out of the crust 
due to transport processes (e.g., see Gourgouliatos and 
Gumming 2015, and references therein). Even if taken at 
face value, it is unclear if these helds can reach the viral 
level of ^ 10^® G in the core. However, as we have argued 
above, this is anticipated by simple astrophysical dynamo 
processes during the formation of the neutron star, fol¬ 
lowed by subsequent cooling in (quasi-)hydrostatic equi¬ 
librium. 

An interesting possibility is that, depending on their 
initial conditions, different neutron stars might have dy¬ 
namos of varying efficiency during their formation pro¬ 
cess. Therefore, future observations of Mass-Radius rela¬ 
tionship of neutron stars may exhibit an intrinsic scatter 
due to variations in internal magnetic held. An even 
more exciting possibility is hnding a surface measure of 
this internal magnetic held (e.g., through modelling the 
braking indices), which could then correlate with (and 
thus effectively reduce) this scatter, following our sim¬ 
ple scalings (Figs. 8-9,13-14). In particular, the sign of 
the correlation will be indicative of the (force-free v.s. 
isotropic) held conhguration. 

In typical MHD turbulence, and in lieu of an active 
dynamo, one expects the magnetic held to relax to a 
force-free conhguration within an Alfven crossing time, 
through subsequent reconnections that conserve helicity 
(Braithwaitc and Spruit 2006; Broderick and Narayan 
2008). Therefore, one may expect the force-free conhg¬ 
uration to be preferred. However, if the interior of neu¬ 
tron stars transitions into type H superconductor, as hy¬ 
pothesized in the literature (e.g., see Haskell et al. 2013, 
and references therein), then it can support a frustrated 
network of flux tubes, which do not reconnect (Brandt 
1995). Such a conhguration can support tangled (statis¬ 
tically) isotropic magnetic helds. 

Finally, one may further speculate that reconnection 
is not completely halted, but is only signihcantly slowed 
(e.g., due to the presence of a mixture of type I and H 
superconducting phases in the interior). Therefore, the 
held conhguration will slowly transition from isotropic to 
force-free, leading to a reduction in the maximum mass. 
Via this mechanism a subset of isolated neutron stars, 
those with masses between those that can be supported 
by force-free and tangled magnetic held conhgurations, 
could undergo a delayed collapse, surviving for a time 
scale set by the evolution of the internal magnetic ge¬ 
ometry. Such an event may present an attractive alter¬ 
nate candidate for energetic, short-time scale phenom¬ 
ena, e.g., short gamma-ray bursts and the recently de¬ 
tected fast radio bursts. 
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APPENDIX 

As mentioned in Section 5, solving the average Einstein equation (20): 

= 8^G( V + (TfiMn), 

along with the energy momentum conservation of the gas: 


requires: 


V T ^ 0 

v gas — '-'5 


= 0 . 


This equation results in the anisotropic TOV equation for the average magnetic field: 

. ^ir). „ , 2 A {B^)/87r 


B 

dr 


= -(eB+Pr)5 + 


where g is given by Equation (25). As is given by Equation (26), and assuming Equation (2), we will have: 


d rPg U-2A 
dr L /3 '' 3 r 


^ 1 
-13 /3^ 3 ’ 


2APg 

j3r 


1 - 2A dPg 2 Pg dA _ A - 2 


3/3 dr 3/3 dr 


3/3 


9 


2APg 

/3r 


( 1 ) 

( 2 ) 

( 3 ) 

( 4 ) 

( 5 ) 

( 6 ) 


and substituting for dPg/dr from Equation (24), and solving for dA/dr, we get: 


dA 

dr 


rU(A- b— 

o ' V 0/70 


12 


2'Pn 


3A 

r 


( 7 ) 
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Figure 15. A{r) for a neutron star of radius i? = 10 km with /3 = 1/3. A remains very close to zero (isotropic field) and more than —1 
until the edge of the star where it drops to lower values. 



R (km) 


Figure 16. Solid curves are the force-free mass-radius relationships already shown in Figure 12. The dashed curves are the results found 
for the self-consistent anisotropy calculation presented in the Appendix. 


which simplifies to: 





gr^dr. 


2'Pgi 


( 8 ) 


The suitable boundary condition for A is determined by the fact that in the anisotropic TOV equation for the 
magnetic field (4), the second term on the right side is well defined only if A —i 0 when r —i 0. Equation (8) solved 
along with the TOV equation for the gas, uniquely determines A(r). As we see in Figure 15, A remains very close 
to zero (isotropic field) and more than —1 until the edge of the star where it drops to lower values. This behaviour 
of A makes at most a 9% change (in case of /3 = 1/3 and smaller changes for larger values of /?) to the mass radius 
relationships previously found as can be seen in Figure 16. It is important to note that values of A < —1 are unphysical, 
implying imaginary magnetic field strengths. This raises some concern regarding the self-consistent solution near the 
stellar surface; for the case shown in Figure 15, roughly 20% of the mass is located in this unphysical regime, implying 
that the mass estimates from the self-consistent solutions are marginal over-estimates. 











